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摘 要 :采用 大 涡 模 拟 方法 计算 Re =2 x10’ 三 维 不 等 直径 串 列 圆柱 (d[D<1) 绕 流 问 题 。 结 果 显 
示 , 处 于 双 涡 脱落 流 态 时 , 随 着 串 列 圆柱 间距 增加 ,上 游 圆柱 量 纲 为 一 的 涡 脱 频率 值 Si 总体 上 升 ， 
而 下 游 圆 柱 量 纲 为 一 的 涡 脱 频率 值 St, 存 在 先 下 降 后 上 升 的 变化 规律 。 在 圆柱 间距 较 小 的 情况 下 ， 
35 随 着 囊 列 圆柱 间距 的 增加 而 减 小 , 量 纲 为 一 的 涡 脱 频率 比值 直径 比 与 间距 比 之 间 近 似 满足 Si,/ 
St, cc (L/D) ~ (d[D) 的 霸 指 数 关系 ;在 圆柱 间距 较 大 的 情况 下 ,圆柱 间 时 均 流 向 速度 提高 并 趋 近 
主流 区 速度 ,Si, 随 间距 比 增加 而 上 升 。 在 较 小 直径 比 囊 列 圆柱 情形 下 ,下 游 圆 柱 量 纲 为 一 的 涡 脱 频 
率 St, 可 下 降 至 更 低 的 临界 拐点 ,从 而 产生 “次 谐 波 涡 脱 锁定 ”现象 。 
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Abstract;The large eddy simulation method is used to simulate the three-dimensional flow past tandem 
cylinders with unequal diameters (d/D<1) at Re 22 x10?. In co-shedding regime , with tandem cylinder 
spacing increasing , dimensionless vortex shedding frequency of upstream cylinder ( St,) keeps rising, and 
the dimensionless vortex shedding frequency of downstream cylinder (St,) has the change of first falling 
and then rising. At small spacing, St, decreases with increased spacing,and dimensionless vortex shedding 
frequency ratio, diameter ratio and spacing ratio are in quantitative power exponent relation , St,/St, oc ( L/ 
D) ""(d/D). M larger spacing, the time-averaged streamwise velocity between cylinders rises and ap- 
proaches the main stream , and St, increases with increased spacing. In tandem cylinders with smaller diam- 


eter ratio, St, moves to a lower inflection point , resulting in *sub-harmonic lock-in” phenomenon. 
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串 列 圆柱 绕 流 问 题 具有 广泛 的 工程 应 用 背 
RU ,例如 航空 飞机 的 起 落架 结构 ,海洋 平台 的 立 管 
和 海底 管线 ,以 及 土木 工程 的 高 层 建 筑 群 等 。 实 际 
工程 中 , 串 列 圆柱 系统 一 般 由 主要 圆柱 和 附属 圆柱 
组 成 ,圆柱 直径 存在 差异 ,而 圆柱 间 的 相互 干扰 作用 
导致 多 圆柱 绕 流 问题 较 单 圆柱 情形 更 加 复杂 ; 因此， 
研究 不 等 直径 圆柱 绕 流 问题 更 具有 普 适 性 和 研究 
价值 。 

针对 串 列 圆柱 的 研究 已 持续 了 几 十 年 ,早期 的 
研究 对 象 一 般 集中 在 等 直径 圆柱 绕 流 问题 上 "1。 
Zdravkovich F 1977 年 首次 给 出 了 实验 结果 并 定义 
其 流 态 类 别 。 文 献 [3-5] 陆 续 研 究 了 雷诺 数 介 于 107 
-G0' 量 级 等 直径 串 列 圆柱 间距 的 影响 。 最 近 几 年 ， 


j 以 及 流体 力 的 影响 。 
CN 近年 来 , 随 着 研究 的 深入 ,不 少 学 者 也 逐步 关注 
不 等 直径 串 列 圆 柱 绕 流 问题 。 其 中 ,在 上 游 圆 柱 直 
E) 相对 于 下 游 圆 柱 直 径 (D) 更 小 的 流动 问题 
(ED <1) 中 ,由 于 流 场 会 呈现 出 丰富 的 流动 状态 ， 
因 奉 得 到 了 较 多 的 关注 。 与 等 直径 串 列 圆柱 绕 流 问 
题 相 比 ,不 等 直径 串 列 圆柱 绕 流 问题 在 流动 特征 上 
姨 和 看 相似 性 也 有 独特 性 。 一 方面 ,不 等 直径 串 列 的 
流 纵 类 别 与 等 直径 相似 ,也 可 分 为 单一 钝 体 绕 流 、 剪 
切 层 再 附着 和 双 涡 脱落 等 流 态 5 。 另 一 方面 , 串 列 
圆柱 直径 比 的 改变 产生 了 不 同 于 等 直径 串 列 圆柱 绕 
流 的 涡 脱落 现象 ,不 等 直径 串 列 的 下 游 圆 柱 具有 两 
个 涡 脱落 频率 值 ,分 别 对 应 上 下 游 圆 柱 各自 的 涡 脱 
现象 ,具有 研究 的 必要 性 。 

在 目前 的 不 等 直径 圆柱 绕 流 问题 研究 中 , 主要 
采用 实验 测量 和 数值 模拟 两 类 方法 。 实 验 研 究 中 ， 
不 等 直径 串 列 圆柱 的 雷诺 数 以 下 游 圆柱 直径 作为 特 
征 长 度 , 介 于 10’ 210 820077, 文献 [ 9] 实验 研 究 
Re =8 x10? ~2.42 x10* 的 不 等 直径 串 列 圆柱 ,发 现 
直径 比 d/D =0.6,0.65 的 串 列 圆 柱 处 于 双 涡 脱落 流 
态 时 ,上 下 游 圆 柱 的 涡 脱 频率 之 间 满 足 两 倍 的 数量 
关系 ,并 将 这 种 频率 关系 定义 为 “次 谐 波 涡 脱 锁定 ” 
现象 。 实 验 研究 中 发 现 这 类 现象 存在 于 特定 直径 比 


的 串 列 圆 柱 绕 流 问题 中 ,但 未 给 出 这 类 现象 与 圆柱 
直径 比 的 关系 ,因而 可 以 通过 增加 圆柱 直径 比 和 间 
距 比 来 探究 “次 谐 波 涡 脱 锁定 ”现象 产生 的 条 件 。 

与 实验 研究 相 比 ,数值 方法 主要 对 二 维 绕 流 问 
题 开展 模拟 , 且 雷 诺 数 一 般 介 于 10” - 107 zn] 05, 
对 雷诺 数 大 于 10 的 研究 较 少 。 例 如 ,文献 [12] 计 
算 了 Re =2 x10’ 不 同 间距 比 时 下 游 圆柱 的 涡 脱 频率 
值 ,文献 [ 13 ] 通过 谱 方法 研究 Re =10° ~ 10? 串 列 贺 
柱 的 临界 间距 值 ,文献 [14] 采 用 二 维 大 涡 模拟 分 析 
Re =10° ,间距 比 L/D =1.75 ~2.5 的 不 等 直径 串 列 
圆柱 的 洲 调 脱落 形态 。 这 些 研 究 都 给 出 了 不 同 间距 
比 下 串 列 圆 柱 的 具体 流 态 ,并 强调 处 于 双 涡 脱落 流 
态 的 下 游 圆 柱 具 有 两 个 涡 脱 频率 值 的 特点 。 然 而 ， 
圆柱 绕 流 的 直接 数值 模拟 结果 表明 5 , 当 Re = 
190 时 , 流 场 已 经 表现 为 三 维 流动 , 且 随 着 雷诺 数 增 
加 ,圆柱 尾 迹 将 发 展 为 不 同 的 二 次 涡 街 。 不 等 直径 
串 列 圆 柱 处 于 双 涡 脱落 流 态 时 ,上 下 游 圆 柱 都 产生 
涡 脱落 现象 ,因而 三 维 计算 结果 才 更 符合 实际 流动 。 
由 于 不 少 学 者 “验证 了 三 维 大 涡 模拟 方法 的 可 行 
性 ,如 文献 [20] 选 择 合适 的 亚 格子 模型 得 到 了 与 实 
验 值 接近 的 流 场 脉动 结果 ,此 外 ,对 于 亚 临 界 雷 诺 数 
的 单 圆 柱 绕 流 , 特 别 是 Re =3 900 大 涡 模 拟 结果 与 直 
接 数 值 模拟 结果 取得 很 好 一 致 性 。 

综 上 所 述 ,为 了 更 准确 地 对 处 于 双 涡 脱落 流 态 
的 不 等 直径 串 列 圆柱 的 频率 特征 展开 研究 ,有 必要 
采用 对 该 问题 的 三 维 模型 进行 变 直 径 比 和 间距 比 的 
数值 模拟 ,并 在 此 基础 上 对 流动 特征 进行 更 深入 的 
分 析 与 讨论 。 因 此 ,本 研究 通过 三 维 大 涡 模 拟 方法 ， 
考察 亚 临 界 雷 诺 数 范围 内 Re =2 x 10 不 等 直径 串 列 
圆柱 绕 流 问题 ,通过 改变 直径 比 d/D =1,0.8,0.6， 
0. 5,0. 4 FEJE L/D =2 ~ 10, 重 点 关注 处 于 双 涡 
脱落 流 态 时 上 下 游 圆柱 涡 脱 频 率 值 随 间 距 比 增加 的 
变化 规律 ,探究 该 雷诺 数 下 的 串 列 圆柱 绕 流 的 “次 谐 
波 涡 脱 锁定 ”现象 ,并 通过 时 均 流 场 和 瞬时 流 场 云图 
分 析 流 动 特征 。 


1 数值 方法 


1.1 ”模型 设置 
建立 不 等 直径 串 列 圆柱 三 维 计算 域 ,借鉴 参考 
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文献 [21] 的 设置 ,如 图 1 所 示 。 上 下 游 圆 柱 的 直径 
分 别 用 d 和 D 表示 ,入口 边界 距离 上 游 圆 柱 中 心 
15D, 横 向 边界 距离 圆柱 中 心 10D ,满足 阻塞 比 小 于 
5% ,圆柱 展 向 长 度 为 4D, 保 证 圆柱 绕 流 的 三 维 流动 
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圆柱 绕 流 的 数值 模拟 工作 ,本 研究 使 用 的 亚 格 
子 模 型 为 Smagorinsky-Lilly 动力 模式 。 采 用 有 限 体 
积 方法 对 上 述 方 程 进行 数值 离散 ,空间 离散 格式 上 ， 
对 流 项 使 用 二 阶 迎风 格式 ,扩散 项 使 用 中 心 差分 格 


充分 发 展 。 圆 心 距 离 用 二 表示 ,为 保证 下 游 圆 柱 尾 
迹 充 分 发 展 , 出 口 边界 距离 下 游 圆 柱 圆心 的 距离 大 
于 20D。 通 过 改变 上 游 圆柱 直径 和 圆心 距离 ,可 以 
改变 不 等 直径 串 列 圆柱 绕 流 问题 的 直径 比 (d/D) 和 
间距 比 (ZXD) 。 对 于 边界 条 件 ,圆柱 表面 为 无 滑 移 
壁面 ,入口 边界 为 速度 进口 条 件 ,出口 边界 为 自由 出 
流 条 件 , 横向 边界 为 对 称 条 件 , 展 向 边界 为 周期 


条 件 。 


图 1 模型 尺寸 和 边界 条 件 


Fig. 1 Model size and boundary conditions 
C, 
1/2; 数值 算法 


民 凡 单个 圆柱 的 特征 雷诺 数 (Re) 来 看 ,本 研究 计 
论 的 不 等 直径 串 列 圆柱 的 Re 范围 为 8 x 10* <Re< 
230 ,此 时 加 柱 尾 迹 已 发 展 为 满 流 ,因而 采用 三 维 
大 渔 模拟 方法 进行 数值 模拟 。 大 涡 模 拟 方法 通过 解 
析 实 尺度 涡 结构 ,使 用 模型 描述 小 尺度 涡 结构 求解 
整 很 流 场 。 大 尺度 涡 结构 由 具体 流动 情况 所 决定 ， 
对 拍 体 流 场 影 响 较 大 ,因而 通过 求解 Navier-Stokes 
方程 得 到 流动 物理 量 值 ;小 尺度 涡 结构 具有 各 向 同 
性 ,适合 通过 建立 模型 计算 。 

三 维 不 可 压 N-S 方程 经 滤波 后 的 控制 方程 为 
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FI p 表示 过 滤 后 的 速度 量 和 压力 量 ; Co 为 来 流速 度 ; 
p 为 流体 密度 ; 为 流体 动力 黏度 ;7; 为 未 知 的 亚 格 
子 尺 度 雷 诺 应 力 ,代表 了 小 尺度 涡 对 大 尺度 涡 的 作 
用 。 通 过 对 封闭 亚 格 子 模型 进行 求解 ,借鉴 前 人 单 


式 , 时 间 离 散 采 用 二 阶 隐 式 格式 。 
1.3 ”模型 验证 


1.3.1 网 格 划 分 

整个 流 场 区 域 网 格 划 分 为 靠近 圆柱 的 区 域 . 过 
渡 区 域 和 远 场 区 域 3 部 分 。 在 两 个 圆柱 周围 生成 
“0” 型 块 结构 ,并 细 化 近 壁 面 边界 层 网 格 ,网 格 增长 
率 为 1. 058 ,保证 垂直 于 壁面 的 第 一 层 网 格 y” <1。 
此 外 ,分 析 了 展 向 网 格 节 点 数 对 圆柱 平均 阻力 系数 
的 影响 ,如 表 1 所 示 。 综 合 考虑 计算 资源 消耗 和 计 
算 结果 后 ,选择 Mesh 2 网 格 设置 并 在 展 向 方向 对 二 
维 网 格 进 行 均匀 拉 伸 ,最 终 形成 三 维 计算 网 格 , 二 维 
网 格 示意 图 如 图 2 所 示 。 


(a) 整体 网 格 


0 


(b) 圆柱 周围 网 格 
图 2 计算 网 格 二 维 示意 图 
Fig.2 Two-dimensional computational mesh 
表 1 不 同 网 格 设置 下 的 串 列 圆柱 绕 流 结果 
(d/D -l,L/D =2) 
Tab.1 Results of flow past tandem cylinders under 
the setting of different grids (d/D 21,L/D 22) 


网 格 设置 。 网 格 1 网 格 2 网 格 3 直接 数值 模拟 
N 15 25 35 
Ca 0.98 | 0.97 0.97 0.99 
Cp 0.13 0.18 . 0.19 0.17 


本 研究 需要 计算 不 等 直径 串 列 圆柱 在 变 间 距 下 
的 流 场 特征 ,在 控制 网 格 总 量 的 情况 下 ,网 格 单元 数 
量 会 随 间 距 比 适量 递增 。 例 如 ,直径 比 d/D =1, 间 
BEEE L/D 22 ~6 串 列 圆柱 的 网 格 单元 数量 从 404 万 
增加 至 603 万 ,具体 的 网 格 单元 数量 如 表 2 所 示 。 


此 外 ,考虑 串 列 圆柱 直径 比 对 流 场 特征 的 影响 ,圆柱 
直径 比 越 小 时 ,网 格 单元 数量 越 大 ,总 体 而 言 ,本 研 
究 计算 网 格 单元 数量 介 于 404 万 和 760 万 之 间 。 
R2 网 格 单元 数量 (d/D =1) 
Tab.2 Meshelement numbers( d/D 21) 


L/D 2 3 4 5 6 


网 格 数量 4040440 4561312 5237104 5828912 6029 776 


1.3.2 算 例 验证 

首先 ,通过 计算 Re = 10 EEH L/D =2,3.5,4 
的 等 直径 串 列 圆柱 绕 流 问题 验证 数值 方法 。 取 圆柱 
涡 脱 落 过 程 稳定 后 的 30 ~ 40 个 洲 涡 及 周期 为 采样 
时 间 ,对 圆柱 升力 系数 时 间 序 列 进行 传 里 叶 变换 得 
到 量 纲 为 一 的 涡 脱 频率 St, 2 f D/U, „St, =f D/ Uo 六 
FEFA E PUE FE WEA ARE , X i A 8s Je 


= 表 3 等 直径 串 列 圆 柱 流 动 参数 对 比 
V Tab.3 The comparison of flow parameters of 
e tandem cylinders with equal diameters 
= [21 [24 究 
E E 
e 大 涡 模拟 。 ”直接 数值 模拟 大 涡 模 拟 
co Re 10? 10? 10? 
CN 
L/D22 á s 0.165 0.172 0.173 
N , 
=u a 0.98 0.99 0.97 
> 
又 Ca 0.15 0.17 0.18 
1/9535 St 0.143 0.154 0.152 
c Ca 0.86 0.85 0.86 
i Ca 0.17 0.16 0.19 
© 
LD=4 S 0. 188 0.181 0.186 
Ca 1.11 1.01 1.08 
Ci 0.53 0.56 0.55 


de 3 对 比 了 本 研究 与 其 他 学 者 的 数值 模拟 结 
果 , 由 于 等 直径 串 列 圆 柱 Sr, = Sc ,表格 中 串 列 圆柱 
的 量 纲 为 一 的 涡 脱 频率 值 都 用 St 表示 ,本 研究 结 
与 前 人 的 大 涡 模 拟 结果 “和 直接 数值 模拟 结果 
基本 吻合 。 定 量 来 看 ,本 研究 的 Sc 值 与 直接 数值 模 
拟 的 结果 相近 ,Si 值 最 大 相对 误差 为 2.7% ,说 明 本 
研究 采用 的 网 格 和 算法 都 是 比较 可 靠 的 。 图 3 为 串 
列 圆柱 在 间距 比 L/D =4 时 基于 量 纲 为 一 的 Q 等 值 
面 (0Q”=1) 的 量 纲 为 一 的 流向 涡 量 ( 玉 “) 云 图。 均 
Fi CU / D) "对 涡 量 进行 量 纲 归 一 化 ,能 够 看 到 圆柱 
间 和 下 游 圆 柱 尾 迹 中 都 具有 复杂 的 展 向 涡 管 结 构 。 


一 
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图 3 直径 比 wD=1, 间 距 比 L/D =4 串 列 
圆柱 流 场 W. KEC =1) 
Fig.3 W,” contour of tandem cylinders with 
diameter ratio d/D =1 and spacing ratio L/D 24(Q* =1) 

1.4 确定 研究 工 况 

根据 文献 [3 ] 的 分 类 , 串 列 圆柱 的 流动 状态 可 
分 为 单一 钝 体 绕 流 、 剪 切 层 再 附着 和 双 涡 脱落 等 流 
So HWA REH ,处 于 亚 临界 雷诺 数 (10 < Re < 
8 x10') 的 等 直径 串 列 圆柱 ,其 双 涡 脱落 流 态 的 临界 
间距 比 (ZAD)e=3.8。 本 研究 主要 研究 不 等 直径 串 列 
圆柱 处 于 双 涡 脱落 流 态 时 的 频率 特征 ,因而 首先 开展 
串 列 圆柱 临界 间距 比 的 研究 ,确定 工 况 的 具体 流 态 。 
由 于 串 列 圆柱 展 向 中 间 截 面 的 流 场 分 布 能 直观 表现 其 
流动 状态 ,下 文 均 选用 中 间 截 面 的 流 场 分 析 流 动情 况 。 

基于 上 述 串 列 圆 柱 物 理 模型 ,通过 三 维 大 涡 模 
拟 方法 计算 Re =2 x 107 不 等 直径 串 列 圆柱 绕 流 问 
题 ,直径 比 d/D-1,0.8,0.6,0.5,0.4, [REB ELD = 
2~10。 直 径 比 d/D -0.8 串 列 圆柱 的 临界 间距 比值 
(ZLD)c 介 于 2 ~3 之 间 , 当 串 列 圆柱 间距 比 L/D =2 
时 , 流 态 为 剪 切 层 再 附着 流动 ,上 游 圆 柱 剪 切 层 被 拉 
长 ,并 在 下 游 圆柱 前 表面 卷 起 附着 ; 当 串 列 圆柱 间距 
比 L/D =3 时 , 流 态 为 双 涡 脱落 流动 ,上 游 圆 柱 在 圆 
柱 间 卷 起 尝 涡 并 脱落 ,脱落 的 洲 涡 撞击 下 游 圆柱 前 
表面 ,对 应 的 上 游 圆 柱 流 涡 脱落 周期 变化 如 图 4 ~ 
图 5 所 示 。 


(a) t-t, (b) =t; +T,/3 


(c) +27,3 
图 4 HE d/D =0.8 ,间距 比 

L/D =2 上 游 圆 柱 流 涡 脱落 周期 变化 图 

Fig.4 Vortex shedding period of upstream cylinder with 
diameter ratio d/D =0. 8 and spacing ratio L/D =2 


(d) = FT, 
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(c) tt; *2T/3 
图 5 直径 比 w[D=0.8, 间 距 比 


(d) = tT, 


L/D =3 上 游 圆 柱 洲 涡 脱落 周期 变化 图 
Fig.5 Vortex shedding period of upstream cylinder with 


diameter ratio d/D =0. 8 and spacing ratio L/D 23 
随 着 串 列 圆柱 直径 比 的 减 小 ,上 游 圆 柱 占 主流 
的 酸 截 面积 变 小 , 双 涡 脱落 流 态 的 临界 间距 比值 也 


逐渐 缩小 。 由 于 本 研究 起 始 间距 比 为 ZLD =2 ,增幅 
XL, 无 法 确定 串 列 圆柱 临界 间距 比 的 精确 值 ,因而 
娩 曲 。 不 同 直径 比 串 列 圆 柱 双 涡 脱 落 流 态 起 始点 如 


ROHR 
Q OS 双 涡 脱 落 流 态 起 始 间距 比 
© Tab.4 Initial spacing ratio of co-shedding regime 
CN 
—d/D 1 0.8 0.6 0.5 0.4 
A 
BD 5 3 3 2 2 
2S P HA 2x m sz Ht 4 
CRUA AA A M A B) 30 E S E 
c 


e 
13; 量 纲 为 一 的 涡 脱 频率 值 与 间距 比 的 关系 


不 等 直径 串 列 圆柱 处 于 双 涡 脱落 流 态 时 , 量 纲 
为 一 的 涡 脱 频率 值 与 间距 比 的 关系 可 分 为 直径 比 
d/D =1 和 直径 比 d/D <1 两 类 。 直 径 比 d/D -1 Œ$ 
列 圆柱 的 量 纲 为 一 的 涡 脱 频率 值 St, = St, , 且 随 着 间 
距 比 增加 而 上 升 ,并 逐渐 接近 上 游 单 圆柱 绕 流 的 量 
纲 为 一 的 涡 脱 频率 值 (St ) ,如 图 6(a) 所 示 。 

直径 比 wD =0.8,0.6,0.5,0.4 的 串 列 圆柱 处 
于 双 涡 脱落 流 态 时 ,上 游 圆 柱 在 间隔 区 域内 形成 一 
对 稳定 脱落 的 涡 , 撞 击 下 游 圆 柱 的 迎风 面 , 因 直 径 比 
减 小 导致 上 游 圆 柱 对 下 游 圆柱 的 频率 主导 作用 减 
弱 , 下 游 圆 柱 自身 的 洲 涡 脱落 频率 值 出 现 ,此 时 下 游 
圆柱 量 纲 为 一 的 涡 脱 频率 值 St, # Sti o 

随 着 间距 比 的 增加 ,S5 总 体 旦 上升 趋势 并 趋 近 
Sto ,而 5 存在 先 下 降 后 上 升 的 变化 过 程 。 图 6(b) 


c— 
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~ 图 6(e) 中 直径 比 wD =0.8,0.6,0.5,0.4 串 列 圆 
柱 的 间距 比较 小 ,5i, 仅 表现 出 随 间 距 比 增加 而 下 降 
的 规律 。 


0.6 


A St, 
0.5 A 
- - -St=0.26 


LID DD 
(a) d/D-1 (b) d/D—-0.8 
0.6 
>si, 
0:5 —n— St, 


- - -81,70.35 


LID 
(c) d/D-0.6 


0.6 


=- St, 
sh 
- - -St=0.52 


9 3 4 5 6 
LID 
(e) d/D-0.4 
到 6 PAREREA — 9 3 A RE E E e e 


Fig.6 Variation of tandem cylinders dimensionless vortex 
shedding frequency with increased spacing ratio 

考虑 到 随 着 串 列 圆柱 间距 持续 增加 ,下 游 圆 柱 
受到 上 游 圆柱 脱落 涡 的 干扰 影响 会 逐渐 减弱 ,可 能 
会 产生 不 同 的 频率 变化 情况 , 因而 选取 直径 比 
d/D =0.4 串 列 圆柱 ,考察 了 较 广 间距 范围 内 圆柱 频 
率 的 变化 情况 ,发 现 它 展现 出 % 变化 的 两 个 阶段 。 
当 间 距 比较 小 (LD =2 ~6) 时 ,S56, 随 间距 比 增加 而 
逐渐 下 降 ; 当 间距 比较 大 (LD =7 ~10) 时 ,Si, 随 间 
距 比 增加 而 逐渐 上 升 , 即 当 串 列 圆 柱 的 间距 超过 茶 
个 临界 值 后 ,Si 的 变化 规律 发 生 改 变 , 如 图 7 所 示 。 

直径 比 4D =0.4, 间 距 比 LD =3,5,7,9 串 列 
圆柱 的 时 均 流向 速度 云图 如 图 8 所 示 。 间 距 比 
L/D -3,5,7 时 ,Sts 处 于 持续 下 降 阶 段 ,上 下 游 圆 柱 
间 存 在 一 段 低速 区 ,连接 上 游 圆 柱 背 风 面 和 下 游 贺 
柱 迎 风 面 , 且 流 速 小 于 横向 主流 区 速度 。 间 距 比 L7 
D =9 时 ,Si, 处 于 上 升 阶段 , 随 着 间距 比 的 增加 ,主流 
区 流速 较 高 的 流体 与 圆柱 间 流 速 较 低 的 流体 混合 ， 


圆柱 间 的 平均 流速 提高 ,接近 横向 主流 区 速度 ,此 时 
下 游 圆 柱 附近 的 时 均 流 向 速度 分 布 也 与 上 游 圆 柱 
相似 。 
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图 7 Ek d/D-0.4 串 列 圆柱 
TOO 量 纲 为 一 的 涡 脱 频率 值 随 间距 比 增加 的 变化 


Fies Variation of dimensionless vortex shedding frequency 


y= vith increased spacing ratio of tandem cylinders with 
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8 直径 比 d/D =0.4 串 列 圆柱 时 均 流向 速度 云 


Fig.8  Time-averaged streamwise velocity field 


of tandem cylinders with diameter ratio d/D 20.4 
直径 比 d/D =0.4, 间 距 比 LD 23,5,7,9 下 游 
圆柱 的 表面 压力 系数 分 布 如 图 9 所 示 。 由 图 可 知 ， 
下 游 圆柱 渍 止 点 的 压力 系数 随 着 间距 比 的 增加 而 提 
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高 。 间 距 比 LD=5,7 时 ,下 游 圆柱 的 渍 止 点 压力 系 
数 相 近 , 稳 定 在 0.3 左右 ;间距 比 L/D =9 时 ,下 游 圆 
柱 的 请 止 点 压力 系数 达 最 大 值 。 

总 之 , 当 圆 柱 间距 离 足够 远 时 ,上 游 圆 柱 脱 落 涡 
的 影响 减弱 ,圆柱 间 的 时 均 流向 速度 提高 并 趋 近 主 


流 区 速度 ,导致 下 游 圆柱 逐 渐 向 单 圆柱 绕 流 形式 发 
展 ,5i, 变 化 规律 改变 为 随 间距 比 增加 而 上 升 。 


1.0 


0 45 90 135 180 


图 9 直径 比 w[D=0.4 下 游 圆 柱 表面 压力 系数 
Fig.9 Pressure coefficient distribution on 


downstream cylinder with diameter ratio d/D 20.4 


2.2. 量 纲 为 一 的 涡 脱 频率 比值 直径 比 与 间距 比 的 
近似 窜 指 数 关系 

直径 比 d/D =0.8,0.6,0.5,0.4 串 列 圆柱 在 间 
距 比较 小 时 ,Si 处 于 下 降 阶段 ,$i,/5i, 随 着 间距 比 的 
增加 而 递减 。 发 现 不 同 直径 比 串 列 圆柱 的 量 纲 为 一 
的 涡 脱 频率 比值 直径 比 和 间距 比 之 间 存 在 近似 过 
指数 的 关系 , 即 St] St, < (ZLAD) (4/D) , AMA, T 
确定 5= -1/4。 各 种 工 况 的 数据 点 和 拟 合 曲线 如 
图 10 所 示 ,大 部 分 数据 点 都 在 拟 合 曲线 附近 ,其 中 
Sb, 值 最 大 相对 误差 为 -5% , 且 拟 合 曲线 的 决定 系数 
ROgepa 1 ,说 明 能 较 好 地 符合 寡 指 数 拟 合 关系 。 
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Fig. 10 Fitting curve of S1,/8t, and (L/D) "^ (d/D) 
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2.3 ”次 谐 波 涡 脱 锁定 现象 


不 等 直径 串 列 圆柱 实验 中 发 现 ' ,特定 直径 比 
的 圆柱 处 于 双 涡 脱落 流 态 时 ,会 发 生 “ 次 谐 波 涡 脱 锁 
定 ” 现 象 ,St =25St,。 当 串 列 圆柱 处 于 双 涡 脱落 流 态 
时 ,圆柱 间距 离 足 够 远 , 下 游 圆 柱 将 受到 上 游 脱 落 涡 
的 撞击 作用 ,从 而 可 以 将 上 游 圆柱 的 尾 迹 视 作 含 扰 
动 的 来 流 ,5t 对 应 扰动 频率 ,5t, 对 应 扰动 作用 下 的 
圆柱 涡 脱 频 率 , 而 St, 则 代表 无 扰动 的 下 游 圆柱 涡 脱 
频率 。 前 人 研究 表明 "扰动 来 流下 的 圆柱 涡 脱 频 
率 值 与 扰动 频率 有 关 , 可 以 通过 St] St, I St,/ Si R 
示 上 游 圆 柱 尾 迹 对 下 游 圆 柱 的 施加 作用 和 下 游 圆 柱 
受到 的 影响 ,进而 根据 间距 比 增加 时 St, / Sto 和 St,/ 
St 的 变化 确定 “次 谐 波 涡 脱 锁定 ” 工 况 。 

ahr 
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kg 11. KP e A TEE St / Sto St / St 随 间距 比 
增加 的 变化 (Re =8 x 107) 

11 Variation ofSt,/ Sto , St,/St; with increased spacing 


gatio of tandem cylinders in experiment (Re =8 x 10?) 
一 图 11 为 串 列 圆柱 实验 中 St / Sc 和 St,/ Bl i] 
距 厂 增加 的 变化 。 由 于 实验 中 间距 比 定义 为 Ld 
( 最上 游 圆柱 直径 为 特征 长 度 ) , 当 串 列 圆柱 直径 比 
较 末 时 ,间距 也 较 小 ,不 同 直 径 比 的 串 列 圆 柱 的 频率 
值 税 表 现 出 某 种 变化 规律 。 例 如 ,直径 比较 小 (d/D = 
0. 25 ) 串 列 圆柱 的 54, 随 着 间距 比 增加 而 持续 下 降 ， 
直径 比较 大 (4d/D =0.7) 串 列 圆柱 Sc, BER] HE E DH] d 
加 经 历 先 下 降 后 上 升 的 变化 , 蓝 色 标记 代表 Sh 上升 
阶段 。“ 次 谐 波 涡 脱 锁定 ”现象 发 生 在 直径 比 4/D = 
0.6 串 列 圆柱 频率 比 与 直线 St /St, = 2 的 交点 位 置 ， 
如 空心 标记 所 示 。 

本 研究 计算 的 串 列 圆柱 St, / St St / St, BG TR] RE 
比 增加 的 变化 如 图 12 所 示 。 按 其 变化 规律 可 以 分 
成 直径 比 d/D 21 MEIH d/D <1 MX., B 12(a) 
中 随 着 间距 比 增加 , 直径 比 d/D =1 串 列 圆 柱 频率 
比值 沿 着 直线 S/S = 1 逐渐 上 升 。 直 径 比 4/D < 
1 的 串 列 圆柱 St,/54 先 降 后 升 的 变化 规律 在 图 
12(a) 中 表现 得 更 加 直观 ,箭头 方向 代表 下 降 阶 段 的 
变化 趋势 , 随 着 串 列 圆柱 直径 比 减 小 ,SA/Sto 临 界 揭 


~ AE T ee dH 
C | naX IV 1 ERBTII 


点 下 降 至 更 低位 置 ,在 St1/St。 较 大 时 才 转 变 为 St 
Sto 上 升 阶 段 。 
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Fig.12 Variation of St,/St, , St,/St, with 


increased spacing ratio of tandem cylinders (Re 22 x 10? ) 

“次 谐 波 涡 脱 锁定 ”现象 出 现在 d/D = 0. 6, 
L/D=5;d/D =0.5,L/D =3;d/D =0.4,L/D =2, 如 
图 12(b) 中 空心 标记 所 示 。“ 次 谐 波 涡 脱 锁定 ”现象 
是 否 存 在 受到 串 列 圆柱 直径 比 的 影响 , 当 串 列 圆柱 
直径 比较 小 (wD =0.4,0.5,0.6) 时 ,在 St,/ Sto F TUE 
阶段 ,频率 比值 与 直线 St, / St, =2 产生 交点 ,发 生 "次 
谐 波 涡 脱 锁定 ”现象 ; 当 串 列 圆柱 直径 比较 大 (dw/D = 
0.8) 时 ,St,/Sto 的 临界 拐点 位 置 上 升 ,St,/Sto 在 下 降 
阶段 不 会 与 直线 St/St, =2 产生 交点 ,因此 不 会 发 
生 “ 次 谐 波 涡 脱 锁 定 ” 现 象 。 

St/Sto 上 升 阶 段 如 图 12 中 的 蓝 色 标 记 所 示 。 
若 串 列 圆柱 相距 无 穷 远 , 上 游 圆 柱 尾 迹 脱落 涡 的 作 
用 极其 微弱 ,此 时 可 视 作 两 个 互 不 影响 的 单 圆 柱 绕 
流 。 例 如 , 当 直 径 比 wD =0.4 串 列 圆柱 相距 无 穷 远 
H$ St] St =2. 5 ,St,/Sto™1 0 


3 结论 


本 研究 利用 大 涡 模 拟 方法 对 Re 22 x 10? 三 维 不 
等 直径 串 列 圆柱 (d/ 忆 大 1) 绕 流 问 题 展开 数值 计算 ， 
通过 改变 直径 比 d/D =1,0. 8,0.6,0.5,0.4 和 间距 
HE L/D =2 ~10 ,着 重 研究 串 列 圆 柱 处 于 双 涡 脱落 流 
态 时 的 频率 特征 。 根 据 数 值 计算 结果 ,主要 得 到 以 
下 结论 。 

1) 随 着 间距 比 增加 ,55 保持 总 体 上 升 趋势 ,So 
存在 先 下 降 后 上 升 的 变化 规律 。 当 间距 比较 小 时 ， 
St 处 于 下 降 阶 段 , 量 纲 为 一 的 涡 脱 频率 比值 直径 
比 和 间距 比 之 间 存 在 近似 窜 指 数 的 关系 , Sts/St x 
(L/D)' (d/D) ,b = -1/4; 当 间距 比较 大 时 , 随 着 圆 
柱 间 时 均 流 速 提高 ,St, 随 间距 比 的 增加 而 上 升 。 

2) 次 谐 波 涡 脱 锁定 现象 与 串 列 圆柱 直径 比 


有 关 , 当 直径 比较 大 (dAD =0.8) 时 ,% 的 临界 拐点 
位 置 上 升 ,频率 比值 与 直线 St] St, =2 没有 交点 ,不 
会 产生 “次 谐 波 涡 胶 锁定 ”现象 ; 当 直 径 比 较 小 
(d/D =0.6,0.5,0.4) 时 ,9 临界 拐点 位 置 下 降 , 频 
率 比 值 与 直线 St] St, =2 产生 交点 。 

后 续 研 究 工 作 可 以 继续 围绕 串 列 圆柱 双 涡 脱落 流 
态 的 频率 特征 展开 ,通过 细 化 圆柱 间距 比 ,来 确定 不 同 
直径 比 串 列 圆 柱 5 临界 拐点 对 应 的 临界 间距 比值 。 
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